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Abstract 

Within an effective field theory framework we study heavy-quark-antiquark systems 
with a typical distance between the heavy quark and the antiquark smaller than 
1/Aqcd- A suitable definition of the potential is given within this framework, while 
non-potential (retardation) effects are taken into account in a systematic way. We 
explore different physical systems. Model- independent results on the short distance 
behavior of the energies of the gluonic excitations between static quarks are obtained. 
Finally, we show how infrared renormalons affecting the static potential get cancelled 
in the effective theory. 
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1 Introduction 



It is an experimental evidence that for heavy-quark bound states the orbital splittings are 
smaller than the quark mass m. This suggests that all the other dynamical scales of these 
systems are smaller than m. Consistently with this fact, the quark velocity v in these systems 
is believed to be a small quantity, v -C 1. Therefore, a non-relativistic (NR) picture holds 
This produces a hierarchy of scales: m ^> mv ^> mv 2 in the system 0. The inverse of 
the soft scale, mv, gives the size of the bound state, while the inverse of the ultrasoft scale, 
mv 2 , gives the typical time scale. However, in QCD another physically relevant scale has to 
be considered, namely the scale at which nonperturbative effects become important, which 
we will generically denote by Aqcd- 

The study of the heavy quarkonium properties was initiated long ago, using potential 
models. A large variety of them exists in the literature and they have been on the whole 
phenomenologically quite successful (see for reviews). However, their connection with 
the QCD parameters is hidden, the scale at which they are defined is not clear, and they 
cannot be systematically improved. In spite of this, great progress has been achieved over 
the years by relating the potentials appearing in these models to some static Wilson loop 
operators 0, ||, |(| [TJ. This formulation is particularly suitable in QCD because it enables a 
direct lattice calculation and/or an analytic calculation in a QCD vacuum model ||. 

However, there still remains the question of to which extent this pure-potential picture 
is accurate || [10], [II], [12], |7[] . In order to answer this question it seems necessary to develop 



a formalism where the error produced by using a Schrodinger equation with a potential 
obtained from QCD (e.g. by lattice simulations [O) instead of doing the computation in 



full QCD (e.g. by NRQCD lattice simulations |14[]) can be made quantitative |15|]. In this 
work we provide the first steps towards this goal. 

In QED, non-potential effects (sometimes noted in the literature as retardation effects) 
arise as 0(a 3 ) corrections to the energy levels (Lamb shift). They are due to (ultrasoft) 
photons with an energy of 0(mv 2 ). Their existence is signalled in QED by an infrared (IR) 
divergence in the 0(l/m 2 ) potential. Hence, potential models for QED are only accurate 
up to 0(a 2 ) corrections. In spite of that, the static potential is well defined at any order 
of perturbation theory and coincides with the Coulomb potential. In perturbative QCD the 
situation may be regarded as similar or very different from QED, depending on the object we 
are interested in. On the one hand, non-potential effects in perturbative quarkonium bound 
states appear as O(af) corrections as in QED, and also due to ultrasoft gluons with energy 
of 0(mv 2 ) (v ~ a s ). On the other hand, the static potential suffers from IR divergences 
in perturbative QCD ||, and hence, in this respect, the situation is completely different 
from QED. In any case, potential models for perturbative QCD are only accurate up to 
0(a 2 ) corrections. The complete 0(a 2 ) corrections to the perturbative static potential 
have now been calculated |16|], and very recently we have presented the leading log (/i- 
dependent) contribution to the 0(a^) correction [jEJ. The next perturbative improvement 
in the determination of the bottom quark mass from the Y(IS') system Jl7|, [18| , as well as the 
expected precision at the Next Linear Collider for top pair production near threshold, are of 
this order . Hence a systematic treatment of the potential and non-potential effects is also 



important in perturbative QCD for physical applications. Here, we shall further elaborate 
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on this point. 

A more complicated issue is how nonperturbative effects in QCD influence the above 
facts. Obviously, it will depend on the relative size of Aqcd with respect to the other 
dynamical scales. To our knowledge, only in the situation where m ^> mv ^> mv 2 ^> Aqcd; 
can statements be made neatly. In this situation, which holds for the lowest-lying states 
of very heavy quarkonium, the leading nonperturbative corrections can be parameterized 
by local condensates P0j, which have a non-potential origin. Our goal here is to discuss 
other situations where neat statements can be made, restricting ourselves to the situation 
mv ^> Aq CD throughout this paper, but allowing for any relative size between Aq CD and 



mv 2 . 



It is our general aim to take advantage of the remarkable progress that has been made 
in effective field theories for heavy quarks during the last years 0, ^T], |TT|, 
p7| , [7], [L2| to address the various issues mentioned above in a solid QCD-based framework. 
Since the hard (m), soft (mv) and ultrasoft (mv 2 ) scales are widely separated, two EFTs 
can be introduced by sequentially integrating out m and mv. Upon integrating out the 
hard scale m, non-relativistic QCD (NRQCD) is obtained f2|. Proceeding one step further 



and integrating out the soft scale mv , potential NRQCD (pNRQCD) is obtained [11]. This 
sequence of EFTs has the advantage that it allows disentangling of perturbative contributions 
from nonperturbative ones to a large extent, which is of the utmost importance in QCD. 

In this work we shall further study pNRQCD, as we define it in detail in section 2. 
This effective field theory contains a singlet and an octet field, which depend on the relative 
(r) and centre-of-mass (R) quark-antiquark coordinates and on time; potential terms, which 
depend on the relative coordinate and the momentum (and spin); ultrasoft gluon fields, which 
depend on the centre-of-mass coordinate and on time. Since pNRQCD has potential terms, 
it embraces potential models. Since in addition it still has ultrasoft gluons as dynamical 
degrees of freedom, it is able to describe non-potential effects. Indeed, the QED version of 
it, namely pNRQED, has been shown to correctly reproduce the non-potential effects that 



arise as 0(a 3 ) corrections to the binding energies of hydrogen-like atoms |28| and positronium 



27J . Moreover, pNRQCD provides a new interpretation of the potentials that appear in the 
Schrodinger equation in terms of a modern effective field theory language. The potentials 
are nothing but r-dependent matching coefficients, which appear after integrating out scales 
of the order of the relative momentum in the bound state (~ mv ~ 1/r) or any remaining 
dynamical scale above mv 2 . Being matching coefficients, the potentials generically depend 
on a subtraction point [i. This is related to the IR divergences found in |§. The subtraction 
point fi may be regarded as the IR cut-off, which is necessary to define the QCD static 
potential to any finite order of perturbation theory beyond two loops. In pNRQCD this 
//-dependence is natural and should not cause any concern. Indeed, below the scale mv 
there are still dynamical gluons, which are incorporated in the pNRQCD Lagrangian. No 
gluons of energy higher than /i are allowed. Hence the dynamics of gluons in pNRQCD is 
cut-off by /i, and this \i dependence will cancel the explicit \i dependence of the potential 
terms when calculating a physical quantity^. Once the interpretation of the potentials as 

1 Some potentials may also have an extra IR scale dependence, inherited from the NRQCD matching coef- 
ficients, which cancels with the UV divergences in pNRQCD arising from quantum mechanical perturbation 
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matching coefficients is taken seriously, one may look for the running of the potentials. We 
have calculated the leading contributions to the running of the singlet, octet and mixed 
matching potentials (to be defined later) in perturbation theory. 

We enforce the equivalence of pNRQCD to NRQCD, and hence to QCD, to the desired 
order in the multipole expansion by requiring the Green functions of both effective theories 
to be equal (matching). The matching between NRQCD and pNRQCD can be done order 
by order in 1/m. This is due to the fact that we are integrating out, in addition to gluons of 
momentum ~ mv , quarks with energies k° of the order of mv, which are always larger than 
the kinetic energy k 2 /2m of the heavy quark for k ~ mv, so that the latter can be expanded 
in the quark propagator (this will be made more precise in sections 2 and 3). If we denote by 
A mp any scale smaller than mv, at each order in 1/m, the calculation in NRQCD contains, in 
general, arbitrary powers of A mp /mv, while in pNRQCD the counting in A mp /mv is explicit 
and we can decide at which order in A mp /mv we want to carry out the matching. This 
is nothing but the multipole expansion. Moreover, if we restrict ourselves to the situation 
mv ^> Aqcd, we can, in addition, do the matching order by order in perturbation theory. 
Then at the leading order in 1/m and in the multipole expansion, the pNRQCD Lagrangian 
only depends on the perturbative singlet and octet potentials. The perturbative singlet 
potential at this order may be obtained from the standard calculation of the static Wilson 
loop. The perturbative octet potential may also be obtained from the static Wilson loop 
with suitable operator insertions at the end-point strings. At the next-to-leading order in 
the multipole expansion, the pNRQCD Lagrangian contains in addition mixed perturbative 
potentials, which involve interactions with gluons of energy smaller than mv. Furthermore, 
at the same order, the perturbative singlet potential is no longer given only by the static 
Wilson loop, nor is the perturbative octet potential. This is due to the fact that the static 
Wilson loop contains contributions from scales below mv, which must be subtracted in order 
to have the perturbative static potential properly defined as a matching coefficient, namely 
as an object that has contributions only from scales of 0(mv). The static potential thus 
defined is the relevant quantity to be used in the Schrodinger equation when the next relevant 
scale is mv 2 . 

In the situation where Aqcd 3> vnv 2 , there is still a dynamical scale above mv 2 . In this 
case, pNRQCD is not yet the suitable EFT at the ultrasoft scale ~ mv 2 ; neither is the 
singlet potential of pNRQCD the suitable potential to be used in the Schrodinger equation. 
In order to obtain the relevant Lagrangian which describes the ultrasoft degrees of freedom 
(~ mv 2 ), and hence the relevant potential to be used in the Schrodinger equation, we still 
have to proceed one step further and integrate out scales ~ Aq CD . Of course, this can no 
longer be carried out perturbatively in a s . The scales larger than mv 2 (but smaller than 
mv) give rise to nonperturbative contributions to the new potential. We will discuss this 
situation in section |[ 

We distribute the paper as follows. In section |2| we review pNRQCD in detail. In section 
[3] we discuss in detail the matching between NRQCD and pNRQCD at leading order in 

theory and not from the US gluons. This particular scale dependence would cancel against a scale of 0(mv), 
the typical relative momentum in the bound state, and not against the US scale of 0(mv 2 ). Since these 
effects appear at 0(l/m 2 ) and we shall deal with O(l/m ) potentials only, we will ignore this possibility in 
the present work. 
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1/m and at next-to- leading order in the multipole expansion. We also present the leading 
contributions to the running of the potentials in perturbation theory. In sections |] and |5] 
we discuss the various physical situations that arise depending on the relative size of Aqcd 
with respect to mv 2 . In section ^| some model- independent results on the short-distance limit 
of the gluonic excitations between static quarks are obtained. Section [7] is devoted to an 
outlook and conclusions. In [Appendix A| we discuss the general matching structure, while 
Appendix B| we show how IR renormalons affecting the static potential get cancelled in 



in 



the effective theory. 



2 Theoretical Framework 
2.1 NRQCD 

Integrating out from QCD the hard scale (~ m) while considering almost on-shell heavy 



quarks and antiquarks produces NRQCD 0. As noticed in |23], |26|, at this stage it coincides 
with the Heavy Quark Effective Theory (HQET). It describes degrees of freedom (quarks and 
gluons) with energy and momentum less than a certain cut-off, which is much smaller than 
m and much larger than any remaining physical scale. The NRQCD Lagrangian is written 
as a power expansion in 1/m. The maximum size of each term may be obtained by assigning 
the next relevant scale (usually mv in NRQCD and Aqcd i n HQET) to any dimensional 
operator. For the purposes of this work it is sufficient to use the NRQCD Lagrangian at the 
lowest order in 1/m: 

£nrqcd = ^° V + X f iD° X - \f; v F^\ (1) 

where ip is the Pauli spinor field that annihilates the fermion and \ is the Pauli spinor field 
that creates the antifermion, iD° = ido — gA°. 

2.2 pNRQCD: the degrees of freedom 



Integrating out the soft scale, mv, in (|l]) produces pNRQCD [|TTJ. The relevant degrees of 



freedom of pNRQCD may depend in general on the nonperturbative features of NRQCD. 
In this work we will assume that there exists a matching scale /i such that mv ^> ji ^> 
mv 2 , Aq CD , where a perturbative picture still holds. 

Strictly speaking pNRQCD has two ultraviolet (UV) cut-offs Ai and A 2 . The former fulfils 
the relation mv 2 , Aqcd Ax <C mv and is the cut-off of the energy of the quarks and of the 
energy and the momentum of the gluons (it corresponds to the \x above), whereas the latter 
fulfils mv <C A2 <C m and is the cut-off of the relative momentum of the quark-ant iquark 
system, p.Q In principle, we have some freedom in choosing the relative importance between 
Ai and A 2 . Although it is not necessary for the calculations of this work, it is convenient to 



2 Notice that although, for simplicity, we usually refer to the matching between NRQCD and pNRQCD 
as integrating out the soft scale, it should be clear from this paragraph that the relative momentum of the 
quarks in pNRQCD is still soft, and hence it has not been integrated out. 
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have in mind the choice A\jm <C Ai, which guarantees that the UV behaviour of the Q-Q 
propagator in pNRQCD is that of the static one when virtual energies are present^. 

If we denote any scale below A 1; i.e. mf 2 ,AqcD, ■••> with A mp , we are in a position 
to enumerate the effective degrees of freedom of pNRQCD. These are: Q-Q states with 
energy of 0(A mp ) and relative momentum not larger than the soft scale; gluons with en- 
ergy and momentum of 0(A mp ). Let us define the centre-of-mass coordinate of the Q-Q 
system R = (x x + x 2 )/2 and the relative coordinate r = xx — x 2 . A Q-Q state can be 
decomposed into a singlet state 5(R, r, t) and an octet state 0(R, r, t), in relation to colour 
gauge transformation with respect to the centre-of-mass coordinate. We notice that in QED 
only the analogous to the singlet appears. The gauge fields are evaluated in R and t, i.e. 
A/_ t = A^(R,t): they do not depend on r. This is due to the fact that, since the typical size 
r is the inverse of the soft scale, gluon fields are multipole expanded with respect to this 
variable. 



2.3 Power counting 

In order to discuss the general structure of the pNRQCD Lagrangian, let us consider in more 
detail the different scales involved into the problem 

1 

771 , p , , A m p . 

r 

The variables r and m will explicitly appear in the pNRQCD Lagrangian, since they cor- 
respond to scales that have been integrated out. Although p and 1/r are of the same size 
in the physical system, we will keep them as independent. This will facilitate the counting 
rules used to build the most general Lagrangian. 

With the above objects, the following small dimensionless quantities will appear 

— , — , A mp r <1. (2) 
m rm 

Note that 1/pr or A mp /p are not allowed since p has to appear in an analytic way in the 
Lagrangian (this scale has not been integrated out). The last inequality in Eq. (|2|) tells us 
that r can be considered to be small with respect to the remaining dynamical lengths in the 
system. As a consequence the gluon fields can be systematically expanded in r (multipole 
expansion) . 

Therefore, the pNRQCD Lagrangian can be written as an expansion in 1/m (from the 
first two inequalities of Eq. (^)), and as an expansion in r (the so-called multipole expansion, 
corresponding to the third inequality of Eq. (g)). As a typical feature of an effective theory, 
the non-analytic behaviour in r is encoded in the matching coefficients. 



2.4 pNRQCD: the Lagrangian 

The most general pNRQCD Lagrangian density that can be constructed with the fields 
introduced above and that is compatible with the symmetries of NRQCD is given, at leading 

3 Of course, this is not so for fixed energy and virtual relative momenta. In this case the UV behaviour is 
dominated by the relative momentum kinetic term. 
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order in 1/m and at 0(r) in the multipole expansion, by: 

£ p nrqcd = Tr | St (^o - V a (r) + . . .) S + Ot (iD - V (r) + . . .) O j 

+#V j4 (r)Tr {O f r •ES + S t r-EO} + g^^Tr {O f r • E O + O f Or • e} 



(3) 



All the gauge fields in Eq. (|]) are evaluated in R and t, in particular F^ ua = F^ a (R, t) 
and zD O = id O - g[A (R, t), O]. 

= 9(T)e~ singlet propagator 



/// / diAj^ 



6*(T) I e *^°-^e ^-272 j octet propagator 

V / ab 



| 



l9 V A J^r-E a 



singlet-octet vertex 



I 

b U a 



octet-octet vertex 



gf abc Coulomb octet-octet vertex 



Figure 1: Propagators and vertices of the pNRQCD Lagrangian (Qj. In perturbative calcula- 
tions the octet propagator is understood without the gluonic string, using instead the Coulomb 
octet-octet vertex. Also the singlet-octet and the octet-octet vertices produce three diagrams 
each in perturbative calculations with the Coulomb gauge: one with a longitudinal gluon line, 
one with a transverse gluon line, and one with both a longitudinal and a transverse gluon 
line. 

In order to interpret the matching potentials V s and V Q properly, let us consider, for 
instance, the leading order (in the multipole expansion) equation of motion of the singlet 
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field 

id S 

where we have included also the kinetic energy, since it contributes at the same order in v as 
V s to the physical observables of the system. This equation is the non-relativistic Schrodinger 
equation. V s is the singlet static potential if the neglected terms in the multipole expansion 
are of non-potential type, i.e. if the next relevant scale of the system is the energy scale 
mv 2 ~ do ~ p 2 /m ~ V s (r). If this is not the case, for instance if the next relevant scale is 
Aqcd; the neglected terms in the multipole expansion still give potential type contributions 
i.e. contributions that do not depend on the energy/state of the system and V s is not the 
static singlet potential. Therefore, we can identify the matching potentials V s and V with 
the singlet and octet heavy Q-Q static potential, respectively, only if no other physical scale 
appears between the soft scale mv and the ultrasoft scale mv 2 , i.e. if Aq CD < mv 2 . See 
section [| for a discussion of this situation. If this is not the case, in order to define a 
potential, scales larger than the ultrasoft one have to be integrated out. This situation will 
be discussed in section |[ 
We define 

a v s (r) . . _ (C A ri \ av ( r ) 



r V 2 J r 

where Cf and Ca are respectively the Casimir of the fundamental and of the adjoint rep- 
resentation of the gauge group. In QCD Cf = 4/3 and Ca = N c = 3. Va and Vb are 
the matching coefficients associated in the Lagrangian (|3]) to the leading corrections in the 
multipole expansion. The coefficients V 8 , V Q , Va and Vb have to be determined by matching 
pNRQCD with NRQCD at a scale \x of the order of or smaller than mv, and larger than the 
next relevant scale. Since we have assumed that [i is larger than Aqcd, the matching can be 
done perturbatively. In order to have the proper free-field normalization in the colour space 
we define 

1 T a 
S = -£=S O = ^=O a , (5) 



where Tp — 1/2. S and O are 1/2 x 1/2 tensors in spin space. The Feynman rules for the 
propagators and vertices defined by the Lagrangian @ are shown in Fig. [I]. 

In principle, terms with higher time derivatives acting on the fields S or O could be 
considered in the pNRQCD Lagrangian. These terms are redundant in the sense that one 
could make them disappear from the Lagrangian and still correctly predict all physical 
observables by re-shuffling the values of the matching coefficients. We have chosen the 
minimal form of the pNRQCD Lagrangian, where higher time derivatives are absent. In 
fact, one can always get rid of these terms by systematically using local field redefinitions 
without changing the physical observables (spectrum) of the theory. 



3 Matching at leading order in \jm 

In this section we will discuss how to perform the matching between NRQCD and pNRQCD. 
In particular we will explicitly calculate the matching at leading order in the 1/m expansion 
and up to 0(r 2 ) in the multipole expansion. 
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The matching is in general done by comparing 2-fermion Green functions (plus external 
gluons with energy of 0(A mp )) in NRQCD and pNRQCD, order by order in 1/m and order 
by order in the multipole expansion. This is rigorously justified with the choice of cut- 
offs made in section 2.2. If the soft scale is in the perturbative region of QCD (i.e. it is 
larger than Aqcd), this can be done in addition order by order in the coupling constant, a s . 
In perturbation theory, matching calculations are usually done, for simplicity, using Green 
functions which are not gauge invariant 0, £23| |2B|] . If the soft scale is not in the perturbative 
region (i.e. it is comparable to or smaller than Aqcd) one can still perform the matching, 
but a nonperturbative evaluation of the Green functions is then required. In this case it is 
important to consider gauge invariant Green functions only. We will try to present formulas 
for both situations. However, since we will be concerned with the situation 1/r > Aqcd, only 
the formulas suitable for a perturbative matching (which typically are not gauge invariant, 
but easier to handle) will be used in the actual calculations. 

When matching NRQCD to pNRQCD it could eventually happen that one does not 
directly end up with the pNRQCD Lagrangian in its minimal form and some terms with 
higher time derivatives are needed. These subtleties do not affect the matching at the order 
we are working in this paper. 

3.1 Interpolating fields 

The matching can be done once the interpolating fields for S and O a have been identified in 
NRQCD. The former need to have the same quantum numbers and the same transformation 
properties as the latter. The correspondence is not one-to-one. Given an interpolating field in 
NRQCD there is an infinite number of combinations of singlet and octet wave-functions with 
ultrasoft fields which have the same quantum numbers and therefore have a non vanishing 
overlap with the NRQCD operator. Fortunately, the operators in pNRQCD can be organized 
according to the counting of the multipole expansion. For instance, for the singlet we have 

X t (x 2 , t)0(x 2 , x i; t)V>( Xl , t) = Z l J 2 {r)S{K, r, t) + Z l ^{r)r r • E a (R, t)O a (R, r, t) + . . . , (6) 

and for the octet 

X t (x 2 ,t)0(x 2 ,R;t)T a 0(R,x 1 ;t)^(x 1 ,t) = Z^^O^K, r, t) 

+ Z^(r)rr-E a (R,t)S(R,r, £) + ..., (7) 

where ^ 

0(y,x;t) =Pexp ji(/j£ds(y-x) ■ A(x - s(x - y), t)} . (8) 

As it is clear from Eqs. (^j) and (|7|) these operators guarantee a leading overlap with the 
singlet and the octet wave-functions respectively. Higher order corrections are suppressed 
in the multipole expansion and will not be considered in the applications of the following 
sections. In particular, higher order operators in the expansions (||) and (^) contribute to 
the singlet and octet matching at order 0(r 3 ) or higher in the multipole expansion. In the 
following we will not go beyond effects of order 0(r 2 ). We notice that in the expansions (§) 
and (0) the momentum operator p does not appear since we are performing the matching 
at the lowest order in the 1/m expansion. 
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The strings attached to the quarks on the NRQCD side contain soft contributions. Their 
role in the calculation will be either to produce corrections which are exponentially sup- 
pressed or to change the normalization factors Z . They will not change the pole of the 
correlator. This is discussed in detail in appendix A. We will show explicitly how this works 
in the perturbative singlet matching of section |3.4.1| . Note that no assumptions are required 
on the large time behaviour of the gauge fields as it is often claimed in the literature, (like 
AJ(T — ► oo) = 0, which is, in general, not true). 

Notice that the matching for the octet ([/]) is not carried out using gauge invariant oper- 
ators. In a perturbative matching this is not problematic, since then, as we will we see, we 
expect V , which corresponds to the octet propagator pole, to be gauge invariant order by 
order in a s . However, one may prefer to work with manifestly gauge invariant quantities. 
This is the case if eventually one wants to take advantage of nonperturbative techniques like 
lattice simulations. A possible solution consists in substituting the T a colour matrix on the 
left side of Eq. (|^) by a local gluonic operator H a (FL, t)T a with the right transformation 
properties. The NRQCD operator has then the form 

X f (x 2 , t)0(x 2 , R, t) H(R, t) 0(R, x 1; t)V(x 1; t). (9) 

To be more specific a possible gauge invariant matching condition for the octet reads 

X + (x 2> f)0(x 2 , R, t) B(R, t) 0(R, x 1; t)V(x x , t) = 4^(r)B a (R, t)O a {K, r, t) + . . . , (10) 

where the dots indicate overlaps with operators of higher order in the multipole expansion. 
Note that, if we would have identified H in Eq. (||) with a chromoelectric field, the matching 



condition equivalent to (TO) would have produced an overlap with the singlet operator, —S 



stronger than with the octet one, E a O a . From the NRQCD point of view this is due to 
the fact that in Eq. (|9]) nothing prevents the dynamical gluons associated with the inserted 
operator to be soft just as it happens for the gluonic strings attached to the quarks. 



3.2 Matching at order r° in the Multipole Expansion 

In this section we identify the matching operators relevant for the matching of V s , Z s , V D 
and Z Q at order (l/m)° in the mass expansion and at order r° in the multipole expansion. 



3.2.1 Singlet Matching 

In order to get the singlet matching potential V s and the singlet normalization factor Z s , we 
choose the following Green function: 

I s EE (O|x t (x 2 )0(x 2 ,X 1 )^(x 1 )^ t ( 2/l )0(l/ 1 ,y 2 )x(z/2)|O). (11) 

In NRQCD we obtain 

I a = 5 3 ( Xl - yi )5 3 (x 2 - y 2 ){W n ), (12) 

where Wu is the rectangular Wilson loop of Fig. |2| with edges x\ = (T/2, r/2), x 2 = 
(T/2,-r/2), y x = {-T/2, r/2) and y 2 = {-T/2, -r/2). The brackets ( ) stand for the 
average over the gauge fields and light quarks. 
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v 



Figure 2: A graphical representation of the static Wilson loop. We adopt the convention that 
time propagates from the left to the right. Therefore, the quarks trajectories are represented 
by horizontal lines and the equal-time endpoint Schwinger strings by shorter vertical lines. 



Equation @ states that the leading overlap of the Green function ( |i~2"D in pNRQCD is 
with the singlet propagator, see also Fig. |^. Indeed, in pNRQCD we obtain at the zeroth 
order in the multipole expansion: 

I. = Z s {r)5\^ - yi )5 3 (x 2 - y 2 )e- iTV ^\ (13) 

(W n ) is in general a complicated function of T, but we just want to single out the soft 
scale. This means that we only need to know it in the limit T — > oo, i.e. as a 1/T expansion, 
in the NRQCD calculation. We define 

Un(W a ) = u (r)+i^ + o(± 

Therefore, comparing Eq. (^) with Eq. (|13|), we get 

a Va (r) 



for T — > oo 



(14) 



V s (r) 
\nZJr) 



-c f 

Ui{r). 



u {r) 



(15) 
(16) 



Note that the matching above does not rely on any perturbative expansion in a s . How- 
ever, since we are concerned with the situation r < 1/Aqcd the matching can be done in 
addition perturbatively, i.e. the quantities on the right-hand side of Eqs. (|T^) and ( |T5| ) can 
be evaluated expanding order by order in a s . At leading order in a s we get 

V s (r) = -C F — or a Vs = a s , 



Zs(r) 



(17) 
(18) 



3.2.2 Octet Matching 

In order to get the octet matching potential V a and the octet normalization factor Z a , we 
choose the following Green function: 



xi+x 2 T\ 
^ ,Xi;-)VW 



x^{Vi)<t> (yi 



yi + y2 



T 



^4^,y 2 ;-|hMlo)- (19) 
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NRQCD pNRQCD 

Figure 3: The matching of V s and Z s at next-to-leading order in the multipole expansion. 
On the left-hand side is the Wilson loop in NRQCD, on the right-hand side are the pNRQCD 
fields. The wavy line represents the ultras oft gluon propagator. 



In NRQCD we have 

if = 5 3 ( Xl - yi )5 3 (x 2 - y 2 )(T a W n T b ), (20) 

where the colour matrices are intended as path ordered insertions on the static Wilson loop 
in the points (R,T/2) and (R, -T/2). 

Eq. (|7|) states that the leading overlap of the Green function (E^) in pNRQCD is with 
the octet propagator. Indeed, in pNRQCD we obtain at the zeroth order in the multipole 
expansion: 

if = Z (r)«5 8 (x! - yi )5 3 (x 2 - y 2 )e- iT ^(0(T/2, -T/2)$), (21) 
where the Schwinger line 



f f T / 2 

<j>(T/2, -T/2) = <f>(T/2, R, -T/2, R) = P exp <^ -ig / dt A (R, t) 

I J -T/2 



is evaluated in the adjoint representation. 

As in the singlet case, we only need to know the 1/T expansion of (T a WnT b ). We define 

L b _S^L r=Mr)+ m + 0(^) for T^oo. (22) 



Therefore, comparing Eq. (P0| ) with Eq. (^l|), we get 



Voir) = (^-C F )^l = v {r), (23) 



2 

lnZ (r) = v x {r). (24) 

Again, the formulas above do not rely on any expansion in a s . If in addition the matching 
can be done perturbatively (r < 1/Aqcd), at leading order in a s we get 

v (r) = ^y~ c *)j- or «^ = «s, (25) 

Z {r) = Tp. (26) 

The matching discussed above is gauge dependent. Still, the octet matching potential 
defined above is expected to be a gauge independent quantity at any finite order in pertur- 
bation theory, since it corresponds to the pole of the propagator, i.e. to an eigenvalue of 
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+ — ® + 



NRQCD pNRQCD 

Figure 4: The matching of V and Z Q at next-to-leading order in the multipole expansion. 
On the left-hand side is the Wilson loop in NRQCD with colour matrices insertions, on the 
right-hand side are the pNRQCD fields. 

the spectrum (just in the same way as the pole mass of the quark propagator is a gauge 
independent quantity at any finite order in perturbation theory [|29[| ). At this respect it is 
easy to see on dimensional grounds (by Fourier transforming T to energy space) that the 
string (0(T/2, — T/2)^) does not give contribution to the potential at any finite order in 
perturbation theory, although it does to Z Q , which turns out to be a gauge dependent object. 

The octet matching can also be obtained in a gauge invariant way by choosing, for 
instance, the following gauge invariant Green function 



(y:,^; -|) b'(^; -f ) # (^±iy 2 ; -() *k>l<>>. < 27 > 

Eq. fllOD states that the leading overlap of the Green function (E7I) in pNRQCD is with the 
octet propagator. In a way similar to the gauge dependent matching, we define for T — > oo 

lln (B(R,r/2)i yD B(R-r/2)) =< , (r)+j ^M +0 Q (2g) 

T (B a (R,T/2)0(T/2,-T/2)^ J B b (R,-T/2)) ° V ; T \T 2 J v ; 



and we get 



V (r) = v 9 \r), (29) 
lnZ Bi0 (r) = vt{r). (30) 



We note that V has to be independent of the operator used as far as it has a leading 
projection on the octet state, while in general the octet normalization factor is not. 



3.3 Matching at order r in the multipole expansion 

At 0(r) there are no additional contributions to the singlet and octet matching potentials 
and to the normalization factors. At this order one finds the leading contributions to the 
mixed potentials Va and Vb- For the purposes of this paper they are only needed to be 
known perturbatively at leading order in a s , 

V A (r) = 1, V B {r) = 1. 
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3.4 Matching at order r 2 in the multipole expansion 

At this order we find the next-to-leading contributions to the singlet and octet potentials 
and to the singlet normalization factor. 

3.4.1 Singlet Matching 

The next-to-leading order correction in the multipole expansion to the singlet matching ( p~3"D 
is given by (see Fig. |3|) 

I s = Z s {r)5\^ - yi )5 3 (x 2 - y 2 )e- lTy ^ (31) 

l-fWj(r) f T/2 dt f ^e-^^-^(r-E a (t)0(t,t')^r-E 6 (t'))') , 

JM C J-T/2 J-T/2 J 



X 



where fields with only temporal arguments are evaluated in the centre-of-mass coordinate 
(we will use this notation in the rest of the work). Comparing Eqs. (|T2|) and ( |31~D (for large 
T), one gets the singlet normalization factor Z s and the singlet matching potential V s at the 
next-to-leading order in the multipole expansion. Eq. ( [31]) does not rely on any perturbative 
expansion in a s . However, since we are considering the case where the matching scale \x is 
much larger than Aqcd the equality (|31*|) can be evaluated order by order in a s . Note that 
this happens even when the correlator on the right-hand side of Eq. fl3~l|) in the physical 
system is dominated by nonperturbative effects. We shall illustrate the use of Eq. fl3"I|) by 
calculating the leading log contribution to the static potential. This will be done in two 
different ways. 

The first way consists in a suitable modification of perturbation theory which is suggested 
by the right-hand side of (|3"TD. a s appears in the second term explicitly in V a and V s (also in 
Va at higher orders) and implicitly in the correlator of gluonic operators. We shall expand 
the implicit dependence, i.e. the gluonic correlator, and keep the remaining dependence, i.e. 
V — V s , unexpanded. At tree level (in dimensional regularization) we have 

d 2 * j f d k 



T F (Et(m,*)$$V)) = ~ifC F C A ^ d J ^hke-^\. (32) 

Using the identity (in Euclidean space) 

d d k k n 2ttI 1 7T T(d + n) 



(27i) d (a + k) v v ' T{d/2) {2n) d sm(n(d + n - 1)) T(d + n - v+ l)Y{y) 
we get in the T — > oo limit (d = 4 + 2e) 



rf+Jl — V 



I s = Z s (r)5 3 ( Xl - yi )5 3 (x 2 - y 2 )e~ iTV ^ 

x ( 1 + % TC F -%iy o - V s ) 3 (- + lE + In (V ° ~ \ s) 2 + constant \ 
7i 3 I e 4tt[i z 1 



+cAr\V - V s f |i + lE + In {Y \J^ 2 + constant! + O (I ) ) . (33) 
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a) 



b) 



i i 
i i 
i i 
i i 
i i 



i i 
i i 
i i 
i i 
i i 



Figure 5: 



Graphs contributing to (W n ) giving rise to Ten* In r (V — V s ) and to 



a 3 lnr (V — V s ) terms. Dashed lines represent Coulomb exchanges. 



where Va has been substituted by its tree level value and V a — V s is also meant to be 
substituted by it. Note that the UV divergences can be re-absorbed by a renormalization 
of V s and Z s . Eq. ( p3|) also contains the leading logarithm dependence on /i of V s and 
Z s . Here we want to check explicitly the cancellation of the hx(V — V s ) dependence, which 
plays the role of an IR cut-off, with the right-hand side of Eq. fll2|) by calculating at the 
same order the relevant graphs of (W n )- These are shown in Fig. |5|. The resummation of 
ladders (longitudinal gluon propagators in the Coulomb gauge) amounts to keep V — V s 
unexpanded. The graphs a) were first discussed in || and calculated in The graphs b) 
(and the symmetric ones) are calculated here for the first time. They involve the end-point 
strings and contribute to the normalization factor only. We get for T — > oo 



a) 5-\n(Wu) 

b) S^MWu) 



24 r 7r 
i C F C\ a 3 



V.) V) + H(Yo - V s fr 2 ) + O (± 



T 6 



7T 



ln((K 



K)V) 



Comparing Eq. (O) with Eq. (R^) we get (V — V s = CACt s /2r at tree level) 



V s (r,fi) 
lnZ s (r,fx) 



-C, 



a Va {r) 



two— loops 



C F C\ a s oq 
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r 7r 



In r/i, 



(Ml(r))two-loops + 



C F C\ al 



7T 



In r/i. 



(34) 
(35) 



(■Uo(r)) two -ioops is known [|T^]. See section | for a discussion. Note that the new contributions 
given in Eq. fl3~4j) and ( p5l) would be zero in QED. The fact that ay B depends on the IR 
behaviour of the theory is, therefore, a distinct feature of QCD. Moreover, we stress that in 
order to match the normalization factor (j35| ) it is necessary to take into account contributions 
coming from the end-point strings, which can be considered irrelevant only at order T°, i.e. 
for the matching of the potential. 

The second way of carrying out the calculation is by straightforward perturbation theory 
in a s , and is standard when doing matching calculations (see [f23|, [OJ 0). If we 

expand in a s and in the ultrasoft scale (V — V 8 ), both expressions flT2"D and (^) are IR- 
divergent (Eq. (|3TD is also UV- divergent). By regulating both expressions in dimensional 
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regularization, the calculation in pNRQCD gives zero (there is no scale) , while the calculation 
of the Wilson loop shows up an explicit dependence on the infrared regulator \x via a typical 
In \xr term. In this situation the matching reduces to enforce the following equality order by 
order in a s (for large T) 

Z s (r)e- lV ° T = (W n ) , (36) 

where the exponential has to be perturbatively expanded in a s . In the right-hand side of Eq. 
(j36|) the first IR divergence appears at order 0(«g) at the two-loop level and it only affects 



the normalization factor. It comes from the diagram a) of Fig. |5| without any Coulomb 
exchange except the ones connected to the transverse gluon and from the diagram b) of Fig. 
|5] with only one extra Coulomb exchange besides the ones attached to the transverse gluon. 
The leading log correction to the normalization factor reads 

Z s (r)=C A (l + SZ s (r)) , 5Z s {r) = ^1^1 hir/z. 

2 7T 

We can see that the scale dependence coincides with the calculation previously performed 
in the effective theory (see Eq. (|35D ). At order 0(a^) (three-loops) we have the following 
equality 

C A (iT8Z s {r)^- - iTbv\ = 5(W n ), 

where the computation on the NRQCD side corresponds to the diagrams a) of Fig. |5| with 
only one extra Coulomb exchange besides the ones attached to the transverse gluon and 5V S 
stands for the leading logarithm correction to V s . We obtain for 5{Wn) (up to terms which 
do not contribute to the potential) 

S(W a ) = &TC A °ZQ (C F + %) Inr M . 
it r 6 V 2 / 

We can now obtain the leading logarithm correction to V s , or equivalently to «y s . As expected 
we get the same result as that one reported in Eq. (j3~4D - 

3.4.2 Octet Matching 

The next-to-leading correction to Eq. (|2TD in the multipole expansion is given by (see Fig. 



If = Z (r)5 3 (^ - Yl )5 3 (x 2 - y 2 )e- iTV °^ ( (0(T/2, -T/2)^> (37) 

.2 r T/2 

C 

-T/2 J-T/2 



-|Wj(r) f T/2 dt f dt'e-^ t - t '^' v °\<P{T/2,tf3v • E a '(i)r ■ E b '(t')^t', -T/2)$) 

N r J-T/2 J-T/2 



We have omitted above a term proportional to (see Fig. |B|) and terms which contain op- 
erators like Tr(r l r J [D l , E J ]00^), which should also be included in the pNRQCD Lagrangian 
at order 0(r 2 ), because they neither contribute to the octet matching potential nor to the 
normalization in perturbation theory. Hence formula ( p7[) is less general than the analogous 
formula (|3l|) for the singlet matching potential and cannot be used beyond perturbation 
theory. 
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Figure 6: Octet self-energy graph. 

We use Eq. fl3~T|) in order to calculate the leading log dependence of the octet matching 
potential. From what we have learnt in the previous section, it is sufficient to calculate the 
leading log dependence of the gluon correlator in the right hand side of Eq. fl3T|). Up to a 
color factor the contribution of Eq. (|37|) is equal, at leading order in a s , to the one calculated 
in the singlet case (see Eq. fl3T|)). Proceeding along the same line as for the singlet case we 
obtain 

V {r,ri = (?±- C F ) ^ = Mr))two-ioo PS + (% - C F ) (38) 
\ 2 / r \ 2 / 12 r it 

Up to now (vo(r))two-ioo P s is unknown. Nevertheless Eq. (|38D shows that the ultrasoft 
corrections to V D are of relative order 0(a^). This makes clear that the widespread believe 
that V Q suffers from IR divergences at relative order 0(aD |50[ U (likely due to some mixing 
up with IR divergences of the normalization factor) is not correct. We do not present the 
calculation of the normalization factor Z 01 since, as already discussed in section |3.2.2| , it is a 
gauge dependent object. 



3.5 One-loop Running 

As a byproduct of the next-to-leading order calculations done above, we have got the leading 
running in the matching scale [i of the matching potentials V s and V Q . For completeness we 
give here also the running of Va and Vb- 



a ) 



b) 



4h 



Figure 7: Singlet-octet vertex at tree level a) and at one-loop b). 

The coefficient Va controls the running of the singlet-octet vertex shown in diagram a) 
of Fig. [F[ Diagram b ) shows the one-loop correction to the vertex. The divergent part of it 
gives the one- loop \i dependence of Va' 

8 a 

V A (r,fi) = l + -CA-lnrfi. (39) 

3 7T 

Analogously the coefficient Vb controls the running of the octet-octet vertex shown in 
diagram a) of Fig. ^. Diagram b) shows the one-loop correction to the vertex and the 
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b) 



0, 



■ 



Figure 8: Octet-octet vertex at tree level a) and at one-loop b) (symmetric graphs are 
understood). 

divergent part of it gives the one- loop \i dependence of Vb'- 



V B (r,fi) = 1 + -C A — \arpL, 

3 7T 



(40) 



which turns out to be equal to Va- In Tab. |I] we summarize the running of the matching 
potentials in the pNRQCD Lagrangian at order 0(r) and of the singlet normalization factor 
at leading order in perturbation theory. 




Table 1: One-loop running of the matching potentials of the pNRQCD Lagrangian (Qj and 
of the singlet normalization factor. 



4 pNRQCD: A QCD £ mv 2 

In the situation Aqcd < mv 2 there is no relevant physical scale between mv and mv 2 and 
the pNRQCD Lagrangian of Eq. @ only describes ultrasoft degrees of freedom. All non- 
perturbative contributions are encoded into non-potential terms, while the static potentials 
(singlet and octet) coincide with the matching potentials V s and V a calculated above, i.e. 
they are purely perturbative. 

In particular, for what concerns the singlet static potential, the result ( |34]) improves 
the knowledge of it up to the three-loop leading logarithm order, «glnr/i. Substituting 
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(uo(r)) two-loops by its actual value, we get 



a Vs (r, fj) = a s (r) \ l + (a x + 27 E /? ) 



Air 



+ 



7T 



7s (4ai^ + 2A) + ( y + 4tI J /5 2 + a 2 



16 7T 2 



12 7T 



(41) 



where a s is in the MS scheme, (3 n are the coefficients of the beta function, ai was calculated 
in Ref. |H] and a 2 in Ref. |16[ (see |16] for notation). One of the most relevant consequences 
of Eq. (O) is that ay s is not a short- distance quantity as a s (in the MS scheme), since it 
depends on \i. Actually, it can be better understood as a matching coefficient. Similar 
conclusions would follow for ay o , although it is not known with the same accuracy. 

The pNRQCD Lagrangian can be systematically written in a 1/m expansion. In this 
work we are just studying the l/m° term. Nevertheless, for the sake of completeness, we 
would like to show how the Lagrangian looks like when 1/m terms are included: 



C 



pNRQCD 



\ 4m m Am 6 



VJr) - 



K (1) K (2) 



+ ... S 



+ O t [iDn- 



4m 



e! + p 4 

m 4m 3 



- V {r) 



m 

VP VP 



m 



+ 



O 



+gV A (r)Tr {O f r • E S + S f r • E O} + #^^-Tr {O f r • E O + O f Or • e} 



4 ^ 



(42) 



where P is the momentum associated to the centre-of-mass coordinate. In Eq. ( f£2|) the 1/m 
corrections to Vb and to pure gluonic operators as well as the higher order terms in the 
multipole expansion are not displayed. Let us note that the kinetic energy is unavoidable 
when computing physical observables, since ido ~ p 2 /m ~ V s (r) ~ mv 2 . The different V}™' 
can depend on r, p and the spin. They are purely perturbative quantities, just in the same 
way as V s>0 is and have been computed with some degree of accuracy in the literature. 

With fl42|) we have finally obtained the relevant Lagrangian for ultrasoft degrees of free- 
dom we looked for. Nevertheless, we would like to comment briefly about the emergence 
of nonperturbative effects in the observables. They will appear as non-potential terms. At 
next-to-leading order in the multipole expansion they will be encoded into two-point gluonic 
correlators which appear in formulas like the following when calculating the energy shift, 



T/2 
T J -T/2 



dt I dt\E?(t)<P(t,tr a fE%t'))e 



-i(T/2-t)h aj .i e - 



■i(t-f)h 0j J e - 



i(t'+T/2)h s 



-T/2 



(43) 



where h s = p 2 /m + V s (r) and h = p 2 /m + V Q (r). Notice that the nonperturbative contri- 
butions get entangled with the spectrum of the singlet and octet Hamiltonians (h s and h Q ) 
in a non-trivial way. More information on the gluonic correlator in ( fI5| ) cannot be obtained 
up to now from first principle QCD, but rely on lattice simulations ||32|| , on QCD vacuum 
models [^, 0] or on sum rules H35j . Still, if we further consider the situation mv 2 3> Aqcd 
the nonperturbative dynamics can be parameterized into local condensates. This situation 
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is that one described by Voloshin and Leutwyler ||20|| . It follows from ( ]I3| ) upon realising 
that h s — h Q ~ mv 2 ^> Aqcd will make the exponentials oscillate wildly unless t ~ t', so that 
a local expansion of the correlator is justified. On the other hand, if we are in the situation 
Aqcd ~ mv 2 , we can not further simplify the functional dependence on the nonperturba- 
tive dynamics and the observables depend on non-local condensates. To our knowledge this 
situation has only been explored in a model dependent context in [36 



5 pNRQCD: mv > Aqcd > mv 2 

In this section we study the situation mv ^> Aqcd ^ mv 2 . We shall restrict ourselves to 
pure gluodynamics except in subsection where we briefly discuss the case of QCD with 



light fermions. Hence, Aqcd here should be understood as a scale of the order of the mass 
of the lightest glueball. 

In this situation, the pNRQCD Lagrangian of Eq. (|3]), which, we remember, has been 
obtained from NRQCD by simply integrating out the soft scale mv, does not only describe 
ultrasoft degrees of freedom (~ mv 2 ) but also degrees of freedom associated with the scale 
Aqcd which are still dynamical. If we are only interested in the lower lying states, it is 
convenient to integrate out the scale Aqcd i n order to obtain a suitable effective theory 
for US degrees of freedom only. We will denote this new effective theory by pNRQCD'. 
Since the nonperturbative dynamics is important in this case, it is more difficult to obtain 
model-independent results. Nevertheless, the multipole expansion still holds when matching 
to pNRQCD'. This point is of utmost importance in order to control how nonperturbative 
effects are actually going to appear. If we work at the lowest order in the multipole expansion, 
while the singlet decouples from the octet and the gluons, the octet is still coupled to gluons. 
Let us call gluelumps the adjoint source in the presence of a gluonic field |37|.P| 



H(R, r, t) = H a (R, t)O a (R, r, t). 

We can figure out two situations for the spectrum of these gluelumps (leaving aside the 
perturbative potential). 

i) The nonperturbative dynamics creates a gap Ah ~ Aqcd in the gluelump spectrum, 
where Ah is the mass of a generic gluelump. 

ii) The nonperturbative dynamics creates at most a gap A H < mv 2 <C Aq C d in the 
gluelump spectrum. 

Situation i) is the generic one from a Wilson renormalization group point of view (assuming 
naturalness of the coefficients) and hence the most likely to occur. However, situation ii) 
cannot be ruled out from first principles. We shall briefly study these two situations in the 
following. 

4 We will see in section ^ that these eigenstates share some relation with the energies of gluonic excitations 
between static quarks in the short-distance limit. 
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5.1 A H ~ Aqcd 

As we have just mentioned the pNRQCD Lagrangian of Eq. (|3]) is not the most suitable 
effective theory in order to describe the lower lying states in this situation. In particular, 
since there is a scale (Aqcd) between the soft scale mv and the ultrasoft scale mv 2 , the 
singlet matching potential V s can no longer be identified with the potential which should 
be used in a Schrodinger equation. However, the latter can be obtained by integrating out 
also the scale Aq CD , i.e. by constructing pNRQCD'. Since at scales smaller than Aq CD 
the nonperturbative effects dominate, it is convenient to have in mind a representation of 
both pNRQCD and pNRQCD' in terms of hadronic fields. Let us identify the relevant 
hadronic degrees of freedom of pNRQCD'. Since gluelumps and glueballs in pNRQCD have 
masses much larger than mv 2 they are integrated out and the only ultrasoft field left is the 
singlet. Let us elaborate on this point using the original fields of pNRQCD. The fact that 
all gluelumps have a gap implies that the octet field develops a gap and must be integrated 
out. Hence the only possible remaining fields are the singlet and ultrasoft gluon fields. If we 
assume that the ultrasoft gluon fields have indeed to be kept, terms as the following should 
be considered, 

SflSrV [cEETiiWW} + c BB Ti{B l BJ}) . (44) 

The term above translates in a hadronic representation to a coupling between quarkonia and 
glueballs. However, for pure gluodynamics the mass of the glueball is believed to be ~ Aqcd, 
a scale we have already integrated out. Hence we arrive at a contradiction. We conclude 
that for the case of pure gluodynamics no ultrasoft gluon fields have to be included in the 
pNRQCD' Lagrangian and hence we demonstrate that in this situation the quarkonium 
dynamics reduces exactly to the motion in a potential. Non- Potential effects do not exist. 
Therefore, the pNRQCD' Lagrangian reads 



£ p nrqcd' = Tr |s f (id - ^ - V s '(r) + . . J S j 



(45) 



where the dots indicate higher-order potentials in the 1/m expansion and the centre-of-mass 
kinetic terms. They are irrelevant here and hence will be neglected. 

Vg is the potential describing systems characterized by mv ^> Aqcd 3> mv 2 . It has to be 
obtained by matching pNRQCD with pNRQCD'. Since Aq C d has been integrated out, the 
matching is nonperturbative and V s ' will contain nonperturbative corrections to the already 

calculated (perturbative) singlet matching potential of pNRQCD, VJr, p) = — Cf— — — ^— , 

r 

see Eq. (|4l|) . These nonperturbative corrections can be computed systematically in the 
matching procedure since the multipole expansion still holds (mv ^> Aqcd)- I n particular 
at the next-to-leading order in the multipole expansion (using the results of section |3.4.1| ) 
we get 

V: = -C F ^^- i ^T F V 2 (r)-^- nte-^^)(E a (t)0(t,O)^E 6 (O))(/i). (46) 
r JM C d — L Jo 

/i is the UV cut-off of pNRQCD (/x > Aq C d)- We notice that Eq. flUf) is close to that one 
obtained by Balitsky in |3S| . Since V Q — V s ~ mv 2 Aq C d in order to perform correctly the 
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matching we have to expand Eq. (JiS) in V Q — V s . In particular, this expansion guarantees 
that only the nonperturbative part (which contains a scale) of the non-local gluon condensate 
on the right-hand side of Eq. ([5]) contributes (in dimensional regularization). Therefore, 
we get 

V; = -C F ?^± - i^rT F Vl(r)-^- [°°dt (l - it(V a - V s ) - \t\V - V s f 
r JS C a — 1 Jo l 2 

+i\t\V - V s f + . . .} (E Q (t)0(t, 0)^E fc (0))(/i) non P ert -. (47) 

Each new term of the expansion is suppressed by an additional factor (Vo — K,)/Aqcd- The 
leading term cancels the 0(r 2 ) renormalon of the perturbative potential (see appendix B), 
while the term proportional to (Vo — V,) 3 cancels the ln/i dependence of ay s (r, /i) (see Eq. 
(PI)). In the following subsection we will explicitly show this cancellation assuming for the 
non-local condensate of Eq. ( |47|) a simple model. 

Finally note that, although the nonperturbative correction can be systematically com- 
puted in the multipole expansion, the leading 0(r 2 ) nonperturbative term could be as im- 
portant as the perturbative potential once the power counting is established, if so, it should 
be kept exact when solving the Schrodinger equation. 

5.1.1 A Model with a Massive Gluon Propagator 

In general, the non-local gluon condensate (E a (t)(p(t, 0)^j j E 6 (0)) should reproduce at high 
energies (let us say at the scale /i of the matching between NRQCD and pNRQCD) the per- 
turbative expression (|32|) while at the scale of the matching between pNRQCD and pNRQCD' 
it is dominated by nonperturbative physics. A crude assumption (which however satisfies 
these requests) is that it can be effectively described by a massive gluon, where the mass 
m g ~ Aqcd plays the role of a hard cut-off at the scale of the matching between pNRQCD 
and pNRQCD'. Then, we have 

where Wk = Jk 2 + m? g . Putting Eq. (^) in Eq. ([47]) and using (in Euclidean space) 

r d d k 1 _ {a 2 )^- n T{n-d/2) 
J (2n) d (a 2 + k 2 Y ~ (4tt)# W ' 

we get, after renormalization in the MS scheme and neglecting a s corrections in Va, 

v; = - CF ^Ml + j!i + 3 (K _ Vtfms _ M ln ^ j . m 

The leading nonperturbative correction in the multipole expansions (which is an expansion 
in mv 2 jmv and mf 2 /AqcD) to the static potential is in this model of order r 2 AQ CD with 
the "right" positive slope. The correction of order (V Q — V s ) 3 cancels the ln/z dependence 
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of the perturbative potential (see Eq. ([34])) so that the expression ((49]) is scale independent 
(as it should be in a theory, like pNRQCD', free from ultrasoft corrections). V!, contains, 

a 4 

therefore, a term like — In rm g which shows that, also in a situation where potential-like 

r 

nonperturbative corrections exist, it would have been wrong to interpret the re-summed 
static Wilson loop (a la Appelquist-Dine-Muzinich) as the singlet static potential. 

Finally we stress that the above results show how, in the physical situation mv ^> 
Aqcd 3> mv 2 , nonperturbative contributions cancel the scale dependence of the matching 
potential V s in pNRQCD, giving rise to a scale independent static potential V!,. The price 
to pay is that now the static potential is sensitive to the nonperturbative physics. Since 
we are still in the region r ~ 1/mv < 1/Aqcd the computed nonperturbative corrections 
show up as a power series in r. In principle all powers are allowed. In the model we have 
used the leading correction goes like r 2 . Since some authors (see for instance [^]) claim the 
existence of a linear potential at short-distances, it is worthwhile to mention that such a kind 
of correction is in principle not excluded even in this model. They may show up at higher 
order in the multipole expansions (i.e. suppressed by powers in mv 2 /mv and mf 2 /AqcD)- 

5.2 A H < mv 2 < A QCD 

If some gluelump fields do not develop a gap, they must be retained as low energy degrees 
of freedom. Then, the hadronic pNRQCD' Lagrangian would contain the light gluelump 
fields in addition to the singlet. In terms of the original fields in pNRQCD, pNRQCD' 
would still contain the octet in addition to the singlet field. No ultrasoft gluon fields would 
have to be included for the same arguments as in the previous subsection, and hence non- 
potential effects would not exist. The singlet and octet matching potentials would receive 
further contributions from the integration of Aqcd in a similar way as the singlet matching 
potential did in the previous subsection. Furthermore, at 0(r) in the multipole expansion, 
we would have mixing of the light gluelumps and the singlet. 

It is worth stressing that this situation does not correspond to standard potential models. 
Since we have mixing of different degrees of freedom (gluelumps) with what would be the 
usual quarkonium (the singlet), it rather corresponds to a coupled channel system in quantum 
mechanics. 

5.3 pNRQCD' with light quarks 

Let us sketch how the effective theory we have called pNRQCD' is modified if light quarks are 
taken into account, i.e. for real QCD. We first recall at this point that the nonperturbative 
dynamics of pure gluodynamics and that of QCD with (more than one) chiral fermions 
has important differences. Whereas the former is believed to have a gap (~ Aqcd) i n the 
spectrum (being the ++ glueball the state with lowest energy), the latter is gapless in 
the limit of exact chiral symmetry since it contains the Goldstone bosons associated to its 
spontaneous breakdown. Still, also for QCD with (more than one) chiral fermions it is 
believed (and experimentally confirmed) that there is a mass gap between the Goldstone 
bosons and the rest of the spectrum that is usually denoted by A x , of the order of the p mass 
(while it is not clear in this case whether the glueball is a stable particle). 
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Within this framework, and assuming that A x ~ Aqcd, we can use A x just in the same 
way as Aq CD was used in pure gluodynamics. The analysis below is parallel to the one 
carried out for the latter case. Since mv 3> Aqcd 3> mv 2 , the multipole expansion is valid 
when matching pNRQCD to pNRQCD'. At lowest order in the multipole expansion, while 
the singlet decouples from octet and gluons (and also light quarks), the octet is still coupled 
to gluons, and these to light quarks. Again, we can figure out two situations on the spectrum 
of the gluelumps: 

%) The nonperturbative dynamics creates a gap Ah ~ Aqcd in the gluelump spectrum; 

ii) The nonperturbative dynamics creates at most a gap Ah < mv 2 <C Aq C d in the 
gluelump spectrum. 

Independently of these two situations, there will exist extra hadronic ultrasoft degrees of 
freedom with respect to the case of pure gluodynamics. These are the Goldstone bosons 
associated to the spontaneous chiral symmetry breaking, namely pions and kaons. Hence, 
non-potential effects will exist. This important feature will also survive in the situation 
mv ~ Aqcd- 

In the case Ah ~ Aqcd, some results of pure gluodynamics remain true. All states with 
energies of O(Aqcd) are integrated out. Now, the relevant hadronic degrees of freedom of 
pNRQCD' are, besides the singlet, the Goldstone bosons. Therefore, the hadronic Lagrangian 
in this case will have additional terms, besides Eq. P5|), taking care of the interaction of 
the singlet with the Goldstone bosons. They enter at 0(r 2 ) in the multipole expansion. In 
terms of the original degrees of freedom (gluons and light quarks), some of them would read 
as follows 

Upon hadronisation these terms will give rise to suitable couplings of Goldstone bosons to 
quarkonium. 

In the case Ah < mv 2 , the hadronic degrees of freedom of pNRQCD' are now, besides 
the singlet and the light gluelumps, the Goldstone bosons (again all states with energies 
of O(Aqcd) are integrated out). Now, the hadronic Lagrangian will have additional terms 
taking care of the interaction of the singlet and/or gluelumps with Goldstone bosons. At 
0(r) in the multipole expansion there are no singlet-singlet- Goldstone bosons terms in the 
Lagrangian, but there can be singlet-gluelump- Goldstone bosons, or gluelump-gluelump- 
Goldstone bosons terms. In terms of the original (gluon and light quark) degrees of freedom 
the pNRQCD' Lagrangian would look like pNRQCD, that is with the singlet and octet fields, 
with additional operators involving the light quarks, like 

^^(OT^^TV* 

which, upon hadronisation would give rise to a singlet-gluelump-Goldstone bosons coupling. 

Let us finally note that, in principle, the actual values of Ah in pure gluodynamics or in 
real QCD are different since the octet can interact with light quarks through the gluons. 
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6 pNRQCD: the static limit 



Although pNRQCD is originally designed to study Q-Q systems of large but finite mass, it 
is interesting to study its static limit. In particular, because, in the short- distance limit, it 
corresponds to the study of the gluelumps defined in the previous section. These, in turn, 
correspond to the gluonic excitations between static quarks in the short-distance limit, for 
which there is abundant nonperturbative data available from lattice simulations. We shall 
restrict ourselves to the case of pure gluodynamics. 

In the static limit there are no space derivatives in the Lagrangian, and hence r and R 
are good quantum numbers. The spectrum then consists of static energies which depend on 
r (translation invariance forbids R dependences) and on the only other scale in this problem, 
Aqcd- In this section we will discuss some general properties of the short-distance behaviour 
(the only accessible to our analysis) of these static energies. These can be straightforwardly 
derived in our formalism. 

In the limit Aqcd "C 1/r, the spectrum of the theory can be read from the Lagrangian 
(0). In particular, the leading order solution corresponds to the zeroth order of the multipole 
expansion. At this order the dynamics of the singlet and octet fields decouple. Hence, 
the gluonic excitations between static quarks in the short-distance limit correspond to the 
gluelumps. Depending on the glue operator H and its symmetries, the gluelump operator 
O a H a describes a specific gluonic excitation between static quarks and its static energy Va- 
in NRQCD (as in pNRQCD) gluonic excitations between static quarks have the same 
symmetries of a diatomic molecule. In the centre-of-mass system these correspond to the 
symmetry group (substituting the parity generator by CP). According to it the mass 
eigenstates are classified in terms of the angular momentum along the quark- ant iquark axes 
(\L Z \ = 0, 1, 2, ... to which one gives the traditional names E, II, A, . . .), CP (even, g, or odd, 
u), and the reflection properties with respect to a plane passing through the quark-antiquark 
axes (even, +, or odd, — ). Only the £ states are not degenerate with respect to the reflection 
symmetry. 

In pNRQCD at lowest order in the multipole expansion, besides the already mentioned 
symmetries, extra symmetries for the gluonic excitation between static quarks appear. The 
glue dynamics no longer involves the relative coordinate r. Therefore, the glue associated 
with a gluonic excitation between static quarks acquires a spherical symmetry. In the centre- 
of-mass system gluonic excitations between static quarks are, therefore, classified according 
to representations of 0(3) x C, which we summarize by L, the angular momentum, CP and 
reflection with respect a plane passing through the quark-antiquark axesQ. Since this group 
is larger than that one of NRQCD, several gluonic excitations between static quarks are 
expected to be approximately degenerate in pNRQCD, i.e. in the short-distance limit r <C 
1/Aq CD . We illustrate this point by building up all operators, H, up to dimension 3 and by 
classifying them according to their quantum numbers in NRQCD and pNRQCD in Tab. |2|. In 
Tab. |2] all the operators are intended evaluated in the centre-of-mass coordinates. is not 
displayed since it corresponds to the singlet state. The prime indicates excited states of the 
same quantum numbers. The chosen operators for the II and A states are not eigenstates of 

5 We could also choose P but we prefer this representation for a better comparison with the NRQCD 
quantum numbers. 
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the reflection operator. This is not important since these states are degenerate with respect to 
this symmetry. Let us consider, for instance, the U g . The corresponding odd and even states 
with respect the reflection symmetry are given by (r x Tr{EO})j_ and (r x Tr{EO})|| which 
are the projections orthogonal and parallel with respect to the reflection plane respectively. 
The operators E and D x B (B and D x E) have the same quantum numbers so by default 
one would expect them to project over the same states but we can go beyond that and see 
that they are related by the equations of motion assuring that they will certainly project 
over the same states. From the results of Tab. |2| the following degeneracies are expected in 
the short- distance limit: 



n 

n, 



9 > 



n' ~ 

a 

n' ~ 



A 3 ; 



(50) 



Similar observations have also been made in \W\. In pNRQCD they emerge in a quite clear 
and straightforward way. Moreover, here we can write explicitly the relevant operators. 



Gluelumps 

Qa H a 


L = 1 


L = 2 




r-E ,r ■ (D x B) 








(r-D)(r-B) 




r x E ,r x (D x B) 




K 




r x ((r • D)B + D(r • B)) 






(r x D) l (r x B) j + (r x D) J '(r x B)' 






(r-D)(r-E) 


K 


r - B , r (D x E) 




n u 


r x B , r x (D x E) 




K 




r x ((r • D)E + D(r • E)) 






(r x D)*(r x E) j + (r x D) J '(r x E) i 



Table 2: Operators H for the S ; II and A gluonic excitations between static quarks in 
pNRQCD up to dimensions 3. The covariant derivative is understood in the adjoint repre- 
sentation. D • B and D • E do not appear, the first because it is identically zero after using 
the Jacobi identity, while the second gives vanishing contributions after using the equations 
of motion. 

So far we have just used the symmetries of pNRQCD at lowest order in the multipole 
expansion. In fact we can go beyond that and predict the shape of the static energies by 
actually calculating the correlators 

(0| J ff(R,r,T/2)i/ t (R , ,r , ,-T/2)|0) ~ 5 3 (R - R')<5 3 (r - r 7 ) e~ iTV M r) (51) 
for large T. At leading order in the multipole expansion we obtain 

V H (r) = V {r) + % - ln(^(T/2)0(T/2, -T/2)^H b (-T/2))) , (52) 
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where the T — > oo limit is understood. The general structure of the gluonic correlator is the 
following (the contribution from the continuum is included in the dots) 

(H a (T/2)(f)(T/2, -T/2)^H\-T/2)) nonpevt - ~ h e" iAffT + ti e" iA « T + . . . . (53) 

Since we are in the static limit, 1/T <C Aqcd ~ Ah < A H < . . ., one can approximate the 
right-hand side of Eq. ( |5*2"D for T — ► oo by just keeping the first exponential of Eq. (|5*3|). 
Then we get at leading order in the multipole expansion 

V H (r) = V (r) + A H . (54) 



Formula ( 54j) states that at leading order in the multipole expansion the short-distance 
behaviour of the static energies for the gluonic excitations between static quarks is described 
by the perturbative octet potential plus a nonperturbative constant. The constant A# 
depends in general on the particular operator H, i.e. on the particular gluonic excitation 
between static quarks. Ah is the same for operators identifying states which are degenerate. 
Notice also that Eq. fl5"3| ) can be systematically improved by calculating higher orders in 
the multipole expansion. In particular, one can look at how the 0(3) xC symmetry is softly 
broken to -Dqo/i i n the short-distance limit. 




Figure 9: Energies for different gluonic excitation between static quarks at distance r from 
the quenched lattice measurements of fcp\l> To — 0.5 fm. The picture is taken from ^SljJ . 

Let us compare our results against the best available lattice data (see Fig. [|and |4l|). Let 
us first note that in our formalism we can trivially disentangle the gluonic excitations between 
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static quarks from the states composed by the singlet plus glueballs. This is no longer so in 
lattice simulations where some care should be taken in order to identify these different states. 
We will assume in the following that the states measured in lattice simulations correspond 
to the gluonic excitations. First of all we observe a tendency of all the static energies except 
Il 5 to go up at short-distances which may be interpreted as an indication that they want to 
follow the octet potential shape as required by (|54]). The IT U — £~ and A u — W u — static 
energies also show a strong tendency to form a degenerate doublet and triplet respectively 
at short- distances as precisely predicted in fl5"0"|). However, something strange is observed for 
the g static energies. According to (^) a degenerate doublet and triplet should be observed. 
We can clearly see the doublet Ti g — £+' but, as mentioned before, the shape of H g is not 
compatible with the octet potential. On the other hand A g — H g miss a state IT^ to complete 
the triplet. We suspect that the plotted U g is a superposition of the real U g in the doublet 
and the 11^ in the triplet. This would explain its peculiar behaviour at small r. We expect 
that future lattice simulations will be able to disentangle these two states and confirm our 
predictions. It would also be nice to have lattice data at shorter distances so that fl54|) can 
be further confirmed. It is interesting to notice that also the hierarchy of the states, as 
displayed in Fig. ^|, is reflected in the dimensionality of the operators of table [2| Lower lying 
states are characterized by lower dimensional operators if we use a minimal bases with no 
electric fields. In particular this allows to understand on simple dimensional grounds the 
highly non-trivial fact that the doublet and triplet static energies lay close together for the 
g states but quite far apart for the u ones. 

The results given above allow us to relate, in the short- distance limit, the behavior of 
the energies for the gluonic excitations between static quarks with the large time behavior 
of some gluonic correlators, in particular with their correlation length. It is particularly 
appealing that we can extract results for the exhaustively studied gauge invariant two-point 
correlator for the gluon field strength tensor: 

(O|^(t)0(t,O)^ j F^(O)|O). (55) 

One can parameterize this correlator as a function of two scalar functions: 

(O|E«(t)0(t,O)^E 6 (O)|O) and (O|B a (t)0(t, 0)^ J B 6 (0)|0) , (56) 

with correlation lengths: Te = 1/A# and Tb = 1/A# respectively. From the results of Ref. 



4Tfl displayed in Fig. |9] we can conclude that 

T E <T B . (57) 

So far lattice simulations of the gauge invariant two-point correlator for the gluon field 
strength tensor have not reached enough precision to confirm this behavior |32f . This would 



be a nice cross-check of lattice simulations. On the other hand, recently, a sum rule calcula- 



tion |35[ found evidence in favor of Eq. (^7|) and also a lattice computation of the gluelump 
masses jH]] seem to confirm Eq. ([57]). It would be highly desirable to have more precise 
lattice data in the short- distance limit in order to test our results more quantitatively. 
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7 Conclusions 



We have presented the matching between NRQCD and pNRQCD at leading order in 1/m 
and at next-to-leading order in the multipole expansion. For the singlet state we have 
shown that the three loop IR divergence of the potential and the two loop IR divergence 
of the normalization factor in NRQCD match with the corresponding UV divergences of 
pNRQCD. This is a non-trivial check which confirms that pNRQCD is the correct EFT 
at ultrasoft scales. We have presented the leading log contribution to both the octet and 
singlet potentials at three loops. Incidentally, the two loop octet potential is not available 
in the literature although it should be easily obtained from the intermediate steps of the 
calculations carried out in JIB| . This would provide an independent test to our proof that 
the octet matching potential is free of IR divergences at two loops. We have worked out the 
running of the singlet, octet, Va and Vb potentials and the singlet normalization factor at 
leading order in perturbation theory. We have also shown how IR renormalons affecting the 
static potential get cancelled in the effective theory. 

The static limit has been taken throughout because the matching calculation can be 
done order by order in 1/m and the static limit is the leading term, and not because we are 
particularly interested in the dynamics of static colour sources. One should bear in mind that 
pNRQCD is designed to deal with actual bound state systems made out of quarks of large 
but finite mass. This is why it is important to separate soft from ultrasoft contributions. 
Our definitions of the potentials are not arbitrary but the suitable ones to achieve this 
goal. Furthermore, the sequence of matching calculations QCD — > NRQCD and NRQCD — > 
pNRQCD guarantees that at any desired order in 1/m and in the multipole expansion, the 
results we are going to obtain (for ultrasoft energies) using pNRQCD are exactly the same as 
the ones we would obtain using QCD. What we gain is that calculations in pNRQCD have 
several important advantages in comparison with direct calculations from QCD: (i) each 
term in the pNRQCD Lagrangian has a size which is simple to estimate, (ii) the calculations 
at leading orders reduce to quantum-mechanical calculations, (iii) non-potential effects can 
be systematically taken into account, (iv) perturbative and nonperturbative contributions 
are disentangled to a large extend. 

We believe that pNRQCD clarifies several long standing issues. In particular the proper 
treatment of the IR divergences in the singlet potential in perturbative QCD. If our primary 
object of concern is the energy between two static quarks, then the IR divergences can be 
regulated by the resummation of Feynman diagrams originally proposed in |J . However, if 
our primary object of concern is to find a static potential which can be used in a Schrodinger 
type equation for actual bound state calculations, this is not the right thing to do. The IR 
divergences can be simply cut-off, but the cut-off dependent potentials have to be used in 
a framework where ultrasoft gluons (which are also cut-off dependent) are properly taken 
into account, so that in the calculation of any physical quantity the cut-off dependence 
cancels. This is what pNRQCD does for us. We also would like to emphasize that the fact 
that nonperturbative (lattice) evaluations of the singlet potential as the energy of two static 
quarks are free from IR divergences does not imply that this is the suitable object to be used 
in a Schrodinger type equation, as we have just argued. 

We have outlined the suitable effective field theories for the ultrasoft degrees of freedom 
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when mv 3> Aqcd ^> rnv 2 , both in pure gluodynamics and in QCD with light fermions, 
i.e. real QCD. In pure gluodynamics, we see that non-potential effects do not exist. Still, 
different possibilities appear depending on whether the gluelumps develop a gap of 0(Aq CD ), 
or of 0(mv 2 ) or smaller. In the first case, only the singlet remains as a dynamical degree of 
freedom at the ultrasoft scale, while the integration of the octet field (and gluons) at next-to- 
leading order in the multipole expansion introduces nonperturbative terms in the potential 
which are organized in powers of a s /rAQCD ^ 1 starting with a quadratic potential. In 
the second case, some light gluelumps should be kept as dynamical degrees of freedom at 
the scale mv 2 . This situation does not correspond to standard potential models since we 
have extra degrees of freedom (gluelumps) besides what it would be the usual quarkonium 
(the singlet). In real QCD, new degrees of freedom, the Goldstone bosons associated to the 
spontaneous break down of the flavour SU(3) symmetry (pions and kaons) have to be added 
at the ultrasoft scale. Therefore, to the features considered for pure gluodynamics one has 
to include the fact that, now, non-potential effects do exist. 

Some safe statements beyond perturbation theory can be (and have been) made for 
pNRQCD in the static limit. In particular we have studied, in the short-distance limit, 
what in lattice QCD are usually called gluonic excitations of the static quark potential. The 
symmetries of pNRQCD imply that these states can be classified according to representations 
of 0(3) x C, which are more restrictive than those corresponding to the symmetry group of 
NRQCD with two static sources, namely D^. This allows us to identify some approximate 
degeneracies between states. We have also provided explicit operators for these states. We 
have predicted the shape of their static energies and related these with the correlation length 
of some gluonic correlators. Note that energies of glue in the presence of static quark- 
antiquark pair are commonly studied as a first step in the Born-Oppenheimer treatment of 
hybrid heavy-quark mesons. 

We would like to stress that pNRQCD provides a solid framework for the study of heavy 
quarkonium systems in a model independent fashion. It may be useful both for higher order 
perturbative calculations and for investigating nonperturbative effects in a systematic way. 
We would like to briefly comment below on a few applications that have not been addressed 
in this paper. 

The presented results may become relevant in accurate, model-independent, determina- 
tions of the bottom or top mass. The former is usually determined either by using sum rules 
|42| , [0| or by direct analysis of the T(15) mass |L7|, [18| , while the latter could be determined 
by the study of top pair production near threshold in the Next Linear Collider ||19|| . From 
the perturbative point of view, the running of the singlet potential is the first step towards 
the full calculation of the leading log correction to the next-to-next-to-leading order results 
available at present for the above observables. On the other hand our analysis provides a 
model-independent framework to estimate and parameterize nonperturbative effects. 

Another situation where pNRQCD could provide useful information is when Aqcd ~ mv. 
Most of the observed charmonium and bottomonium states correspond to this situation, and 
hence the results presented here cannot be directly applied. This situation requires that the 
matching between NRQCD and pNRQCD is carried out nonperturbatively. Nevertheless, 
many of the features observed in this work survive in a nonperturbative analysis. The 
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potential terms can still be obtained in an expansion in 1/m, while the multipole expansion 
can also be applied for the ultrasoft degrees of freedom. Moreover, in real QCD the Goldstone 
bosons (pions and kaons) remain dynamical at the ultrasoft scale producing non-potential 
effects. In fact, some of the formulas presented here also hold in that situation. Work in this 
direction is in progress 



Note Added. After completion of this work Ref. |44| appeared, where pNRQCD is applied 



to the calculation of the US contributions to the heavy quarkonium spectrum and to the Q-Q 
production near threshold at next-to-next-to-next-to-leading order. We have also recently 
carried out the complete leading-log next-to-next-to-next-to-leading-order calculation (i.e. 
the order ma^ lna s , in the situation Aqcd <C mal) of the heavy quarkonium spectrum p5f . 
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Appendix A Matching: general remarks 



In order to clarify the matching procedure adopted in the main part of the paper, we will 
study in this appendix the general behaviour of the static Wilson loop, 



-ig i dz^AJz) 
W n = Pe J rxT ^ (58) 

in QED and QCD. In Eq. fl58|) P is the path ordering operator (necessary only for non- 
commuting gauge fields) and the integral is extended on a rectangular path with spatial 
dimension r and temporal dimension T. A graphical representation is given in Fig. 0. 

A.l QED 

The Coulomb gauge is implicitly assumed in the discussion below, although, since we are 
working with gauge invariant quantities, the results remain true for any gauge. 

Let us consider a world with static fermions and soft photons. The static Wilson loop 
can be computed exactly. If we neglect the end-point strings (this is legitimate since, in 
the Coulomb gauge, the A component decouples from the transverse photons) the Wilson 
loop can be calculated giving rise exactly to the Coulomb potential V cou i = —a. jr. Up to 
renormalization, the role played by the the soft end-point strings is to create states with 
a definite number of transverse soft photons. The energy difference between states with 
different number of transverse photons is of 0(mv) by definition (we are only taking into 
account soft photons). These contributions are exponentially suppressed in the pNRQED 
power counting, 0(e~^ v ) (since we take V cou \ ~ mv 2 , -Ephoton = k ~ mv and 1/T ~ mv 2 ), 
with respect to the ground state contribution (no soft photons in the initial/final state). In 
fact, they do not give contributions to the matching at all. This is easily illustrated if we 
consider the Fourier transform of the Wilson loop with respect to T (taking into account the 
9(T) due to the heavy quark propagator). We get 

/ dT9(T)e iET W(V ) * ; + / d ' k v lZ> + ■ ■ ■ ( 59 ) 

J E — Vcoui + J (2ix) 6 E — V COVl \ — k + xe 

Since we are only interested in E ~ mv 2 we are near the first pole and all the other contri- 
butions can be expanded and do not contribute neither to the potential nor to the normal- 
ization. That is, they do not contribute to the matching. Specifically, we want to match the 
expression above to the one obtained in pNRQED (Vqed is the potential in pNRQED) 



E - Vqed + ^ 



In this way one trivially gets Vqed = V cou \. 

Let us note again that transverse (dynamical) photons do not give any contribution to 
the static potential and that states with transverse soft photons are exponentially suppressed 
in the Wilson loop. Vqed is 0(v ~ a) suppressed with respect to its natural, soft, size due 
to the a factor coming from the coupling. 
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A.2 Perturbative QCD 



Most of what has been said above for QED also holds for perturbative QCD but instead of 
being exact it turns out to be only approximate (up to 0(a B ) corrections). 

We again consider a Wilson loop with soft gluons and static quarks in the Coulomb 
gauge. We cannot solve the spectrum exactly now. There are interactions among the longi- 
tudinal and the transverse degrees of freedom. This means that, in this case, and unlike in 
QED, states with a different number of soft transverse gluons interact with each other. The 
analogous of Eq. reads 

s m ^ ETw ^- E -c F :^[ ...) + , + /w g -. + 'ok) + « + - (eo) 

We see that the dependence of ay s on a s (in the MS scheme) is not trivial, a>y s — a s (r) + . . . 
Moreover, a non-trivial normalization factor, Z s (r), appears in QCD. The situation in QED 
with light fermions is similar. 

Although we have chosen a straight end-point string in connecting the static quarks 
another configuration could be used. This change of initial/final states produces a different 
projection on the eigenstates of this world of static quarks and soft gluons, but obviously 
does not change the spectrum. As far as we restrict ourself to the matching and to end-point 
configurations with the same quantum numbers as the ground state the change will only 
affect Z s . 

Finally, let us note that the IR cut-off we have introduced in order to ensure that only 
soft degrees of freedom are considered will eventually appear. In QCD, as we will see, this 
already happens for the static Wilson loop. In QED, due to the decoupling of the transverse 
degrees of freedom, the static Wilson loop does not have IR problems (although they will 
appear at higher orders in the 1/m expansion). 

Summarizing, the role played by the soft scale has been to produce potential terms and 
a normalization factor. It should also be noticed that different end-point strings will only 
change Z s but not the potential. 



Appendix B Renormalons 

In this appendix we study the renormalon ambiguities affecting the matching coefficient V s . 
For our purposes it is enough to introduce the running of a s at one loop within the integral 
in momentum space (which basically means to sum up the leading log of all the bubbles 
diagrams, see Fig. |H]), i.e. 

1 + Po— : m — o 

An /r 

If dimensional regularization is used when matching NRQCD with pNRQCD, the integral 
above runs over all values and Eq. (|61|) suffers from IR renormalons. Let us now see the 
structure of these singularities. It is sufficient to work with a hard UV cut-off such that 
1/r ^> \l ^> Aqcd and concentrate on scales below this cut-off /x so that the multipole 
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Figure 10: Bubble diagrams contributing to the renormalon of the static potential. 



expansion can be applied [46]. We obtain the following general expansion for the potential 
with hard cut-off 

V? = C + C 2 r 2 + ... 



where 



and 



1 r 

C 2 = 2nC F - J 



v- d 3 k 1 
(27r) 3 k2 

1 w d 3 k 



aJk 2 



{2ti) A 3 7T ^ 



n=0 



' P at a (fj,y 

47T 



The coefficients and Dffi are given by 



i / 1 
da; ( In — 

x z 



D m 



dx x 2 fin — 

x 2 



3 n+l 



ni 



(62) 
(63) 

(64) 



(65) 



Due to the non-alternating sign of these coefficients, Co and C 2 are not Borel summable and 
they have an ambiguity of O(Aqcd) and 0(Aq CD ) respectively in their definitions J47J, 



5C ~ Aqcd , SC 2 ~ Aq CD 



(66) 



B.l Leading order in the multipole expansion 

The O(Aqcd) ambiguity (of IR origin) of V s cancels against the O(Aqcd) ambiguity (also of 
IR origin) of the mass self-energy (see Fig. [11]) when computed in the pole scheme (which 
has been implicitly assumed throughout this work) |}48| . Therefore, physical observables are 
free of this renormalon ambiguity. 

In fact when matching NRQCD with pNRQCD also the heavy quark and antiquark 
self-energies appear. One obtains 

2E(. • p = 0) = -iS,C F J ^ d <?> = 4, Cf / |^*<k>) , (67) 
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Figure 11: Bubble diagrams contributing to the renormalon of the pole mass. 



which is zero in dimensional regularization while from a renormalon point of view it has an 
UV and an IR renormalon that cancel each other. Nevertheless, in this respect it is nice to 
see that the full result one obtains when matching NRQCD with pNRQCD, i.e. 

2Z(v ■ p = 0) + V s (r) (68) 

has an UV renormalon while the IR renormalon of each term cancel each other (at lowest 
order in the multipole expansion). This suits the EFT picture where the IR renormalon of 
twice the pole mass, 0(m), understood as a matching coefficient, cancels against the UV 
renormalon of the next order 0(m°) in the 1/m expansion. 

B.2 Next-to-Leading order in the multipole expansion 

The 0(AQ CD r 2 ) ambiguity of V s will be absorbed in the effective field theory. It is expected 
to cancel with the renormalon ambiguity of some nonperturbative effects. Therefore, it is 
different in nature with respect the previous discussed renormalon. The explicit form of 
these nonperturbative effects will depend on the relative size of Aqcd with respect the other 
dynamical scales of the problem (mv 2 , ...) and could be of a potential (r 2 nonperturbative 
potential if the next relevant scale is Aqcd, see section ||) or non-potential nature (Voloshin- 
Leutwyler corrections, ... if Aqcd < mv 2 , see section |j). This is not our problem here 
but rather to check that the IR renormalon of the matching coefficient cancels against the 
appropriate UV renormalon in pNRQCD (see Fig. |T2|), if pNRQCD has to be a sensible 
effective field theory. This is, indeed, the case. 




Figure 12: Bubble diagrams contributing in pNRQCD to the renormalon of the singlet 
self- energy at the next-to-leading order in the multipole expansion. 

The calculation of the UV renormalon of pNRQCD at next-to-leading order in the mul- 
tipole expansion effectively reduces to the term in Eq. ( |3~1~D linear in T. The expression from 
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pNRQCD that should be added to the potential reads 



-^JfV^j-^I rfte-^^)(E a (t)0(t,O)^E fc (O))(/i). (69) 

Working at the lowest non-trivial order in a s and introducing an IR hard cut-off analogously 
to the procedure followed for the potential (for our purposes here it does not matter that 
the IR behavior is no accurately described), we obtain for Eq. fl69|) in the UV regime 

r 2 f d d ~ x k 2 r 2 Cpajji) 3 ~ - (2) (foaMY (7n , 

where the coefficients are given by 

r°° / 1 \ n 2 n 

D^ = l dx X d - 2 (ln-) (71) 

The non-alternating sign of these coefficients makes Eq. ([TOD ambiguous by an amount 
0(Aq CD ), corresponding to the UV renormalon. By closer inspection we can see, as expected, 
that Eq. ( |70"D cancels against the term C^r 2 in (|62|), consequently checking the cancellation 
of the renormalon ambiguities up to the considered order in the multipole expansion. 
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